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The Fermi velocity [vf) associated with the spin-orbit coupling is two orders of magnitude smaller 
for spintronic semiconductors than it is for topological insulators. Both families can be treated with 
the same Hamiltonian which contains a relativistic (Dirac) linear in momentum term proportional 
to vf and a non-relativistic quadratic contribution with Schrodinger mass (m). We find that the 
AC dynamic longitudinal and transverse (Hall) magneto-conductivities are strongly dependent on 
the size oi vf- When the Dirac fermi velocity is small, the absorption background provided by the 
interband optical transitions is finite only over a very limited range of photon energies as compared 
with topological insulators. Its onset depends on the value of the chemical potential (/x) and on the 
magnetic field (B), as does its upper cut off. Within this limited range its magnitude is however 
constant and has the same magnitude of e^7r/(8/i) as is found in topological insulators and also 
in graphene noting a difference in degeneracy factor. The total optical spectral weight under the 
universal interband background is e^F/(8>h)Amv‘f^. In contrast to the known result for graphene 
no strict conservation law applies to the spectral weight transfers between inter and intra band 
transition brought about by variations in the magnitude of the chemical potential when a non- 
relativistic contribution is present in the Hamiltonian whatever size it may have. 

BAGS numbers: 78.20.Ls,71.70.Di,73.25.+i 


I. INTRODUCTION 

Since the initial theoretical discussion and experimen¬ 
tal discovery of topological insulators^”— their topologi¬ 
cally protected metallic surface states have been exten¬ 
sively studied, both because of the new physics involved 
and because of novel functionality that could find appli¬ 
cation in a next generation of electronic devices. Another 
class of materials which is also of great importance for 
possible device applications are those used in studies of 
spintronics^^”— Both classes involve the spin orbit interac¬ 
tion. For topological insulators this interaction is strong 
and the linear in momentum relativistic (Dirac) part of 
the single particle Hamiltonian proportional to the Fermi 
velocity {vp) dominates over a smaller non-relativistic 
(Schrodinger) quadratic in momentum piece character¬ 
ized by a mass (m). By contrast, in the materials of in¬ 
terest for spintronics, the Schrodinger contribution dom¬ 
inates over the weak Dirac contribution. In a first ap¬ 
proximation, the same minimal model Hamiltonian can 
be used to describe both cases, but the magnitude of the 
parameters involved in the Dirac and Schrodinger contri¬ 
bution are very different. 

In presently studied materials, the Schrodinger mass 
can vary from the order of the bare electron mass (rUe) 
to one tenth its value and even less. While the Fermi 
velocity vp is of order 5 x 10^ m/s in topological insu¬ 
lators, in spintronic materials it is less than a few times 
10^ m/s and often much smallerThis large difference in 
the magnitude of the spin orbit coupling between these 
two classes of materials can lead to profound differences 
in their physical properties. For example in topological 
insulators the DC quantum Hall effect shows a quanti¬ 


zation n = 1/2, 3/2, 5/2... while in semiconductors it is 
n = 0,l,2,3..^'iS 

The AC magneto response of materials has been 
widely studied and provides valuable information on 
electron dynamics. Recent related studies include 
graphene^”— silicene^i^^ topological insulators^”— 
Mo5'2^and Weyl semimetalsIn topological insulators 
one is dealing with real electron spin and spin momen¬ 
tum locking has been observed 4i”— For graphene and 
other two dimensional membranes such as silicene which 
displays buckling and planar M 0 S 2 , it is the pseudo spin 
associated with the two sublattices of the honeycomb lat¬ 
tice which is involved. Weyl or Dirac semimetals are a 
three dimensional version of graphene. There have also 
been studies of Kerr and Faraday effects in thin films 
with the breaking of time-reversal symmetr y^^i^^ . Mag¬ 
neto transport studie o^^i^^ with emphasis on the tran¬ 
sition between ordinary and topological insulator states 
described by a massive Dirac fermion model with change 
in sign of the mass term. Studies are also availale on 
magneto-optics of bilayer—, multilayer— graphene and 
graphene on polar substrates^. Finally we mention 
magneto-phonon resonances in graphene studied by Ra¬ 
man spectroscopy^. 

In this paper we study the magneto-optical response 
of a two dimensional electron gas with Hamiltonian 
{H) consisting of a combination of relativistic and non- 
relativistic piece. We feature prominently the effect of a 
small spin orbit coupling when the non-relativistic part 
of H is dominant. We also compare with results, some 
known but many new, that apply in the opposite limit 
of a topological insulator for which the spin orbit cou¬ 
pling dominates. Both dynamic AC transverse {axy{uj)) 


and longitudinal conductivity are considered. 

In section II we provide the formal expressions for the 
magneto-optical conductivity based on our chosen Hamil¬ 
tonian. We also give simplified but approximate formulas 
for the optical spectral weight of the various absorption 
lines which correspond to the underlying Landau level 
(LL) structure created by the magnetic field in the limit 
where the Schrodinger part is dominant. We work to 
leading order of Dirac (Ei) to Schrodinger (Eq) mag¬ 
netic energy scale assuming their ratio to be small. We 
explicitly consider the case when the chemical potential /i 
is much larger than Eq . In section III we provide numer¬ 
ical results for both Re(Jxx{^) and Imaxy{oj) when Ei is 
increased and is no longer much smaller than Eq . Partic¬ 
ular attention is paid to the emergence of the interband 
background associated with optical transition between 
the split helical bands caused by the spin-orbit coupling. 
This background is found to be of constant magnitude in¬ 
dependent of the size of vp but is non-zero only in a very 
limited photon window. We also present results when the 
magnetic field is zero. These greatly help in the physi¬ 
cal understanding of the finite B case. Readjustments 
of optical spectral weight due to variations in chemical 
potential are described. In section IV we turn to the 
limit when it is the Dirac energy which dominates. We 
provide simplified but analytic expressions which show 
the first non-zero corrections to the pure relativistic case 
for the optical spectral weight when a small subdominant 
Schrodinger contribution is also included i.e. Eq/Ei 1. 
Numerical results are provided when our simplified ex¬ 
pression is no longer valid and comparison with the limit 
EijEo 1 is made. In section V we give more details on 
spectral weight redistribution with variations in chemi¬ 
cal potential ji. A summary and conclusions are found 
in section VI. 


II. FORMALISM 

The Hamiltonian on which this work is based takes on 
the form 

Ho = ^— + hVF{kx(Ty - kyCJx) ( 1 ) 

where k is momentum, m is the Schrodinger mass, vp the 
Dirac Fermi velocity and (jy spin Pauli matrices. For 
spintronic materials the first term in (1) is dominant with 
the second small in comparison; while for topological in¬ 
sulator it is the opposite. This is shown schematically 
in Fig.l. The dashed black curve (left frame) is the per¬ 
fect cone which would apply for example to the Dirac 
fermions of graphene and is included here for compari¬ 
son with the heavy solid red curve which illustrates the 
dispersion curve found in topological insulators. In this 
case the subdominant Schrodinger piece in the Hamilto¬ 
nian (1) reshapes the cone into an hourglass figure with 
upper conductance band (red) narrowing in cross-section 
as compared with the Dirac cone and lower valence band 



FIG. 1. (Color online) Schematics of dispersion curves used 
to describe a topological insulator (left hand figure) and the 
opposite limit of weak spin orbit coupling (right hand figure). 

(red) fanning out from the cone giving it a larger cross- 
section. The right hand frame illustrates the other limit 
of dominant Schrodinger quadratic band (black dashed 
curve) with subdominant spin orbit coupling. This leads 
to a splitting of the dashed curve into two bands with 
heavy solid red {Ep ) contained inside the solid blue {E- 
) dispersion curve. 

When a magnetic field B is applied perpendicular 
to the plane of the two dimensional Hamiltonian (1), 
Landau levels (LL) form and the dynamic magneto¬ 
conductivity is given by2& 


iV jiV 5O5O 

^^^ {N,s\p\N',s'){N',s'\MN,s) 

id — En^s + Epp^s' + V ( 2 ^) 

where we have included a small phenomenological con¬ 
stant residual broadening of l/(2r). In Eq. (2) Ip = 
1 / ^/e\B\Jh is the magnetic length with e the electron 
charge, / the Fermi distribution, En,s the Landau level 
energies with eigenfunction |A, s). The matrix elements 
{N, s\ja\N', s') with ja the a' th component of the cur¬ 
rent operator carry the information on the optical selec¬ 
tion rules. The Fermi function /tv,^ = -h 1] 

with P the inverse temperature and ji the chemical po¬ 
tential. Details can be found in reference (26) with their 
Eq. (10) corrected to include an overall additional minus 
sign in the second line. The eigen energies are given by 

Em,s = NEo + (So/2)2 + 2NEf (3) 

for V = 1, 2, 3..., s = E and 

En=o = Eo/2. 


( 4 ) 
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FIG. 2. (Color online) Schematic of possible optical transi¬ 
tions in the limit of small spin orbit coupling. The intraband 
transitions (++ solid black arrows) and (-solid red ar¬ 

rows) have optical spectral weight of order one and remain 
finite when vf — The interband transitions (solid green 
arrows) have weight of order for (+—) and for (—h). 


Here Ei = hvFy/e\B\/h is the Dirac magnetic en¬ 
ergy scale and Eq = is the corresponding 

Schrodinger magnetic energy scale. We begin with the 
limit when, for a given value of the magnetic field 
El < Eq. Reference (26) was exclusively concerned with 
the opposite limit Ei > Eq which applies to topological 
insulators. When appropriate we provide comparisons 
between these two cases. In a later section we also present 
additional new results in the TI limit as well. When 
El 0, the energies En^s {N P s/2)Eo + ^^8NP 
with P = {Ei/Eo)‘^ for N = 1,2,3... and En=o = Eq/2. 
Clearly in this case we are dealing with two separate 
Schrodinger Landau level series with the N = 0 term 
included with En,-\- {N ^ 1). These two series can be 
thought of as originating from the spin degeneracy of a 
simple quadratic band split by a small spin-orbit cou¬ 
pling term which shifts one band up and the other down 
in energy. This is illustrated in Fig. 2 where we show 
the series as red dashed horizontal lines and the 

Eq plus En,-\- as solid black horizintal lines. 


The final form of the equations for the absorptive part 
of the dynamic conductivity Recjxxi^) and Imaxy{uj) re¬ 
spectively for longitudinal and Hall conductivity given 
in Eq. (22) of reference (26) will be our starting point. 


They are 


F{0,s) 


En 


X [5{huj - EoH{0, s)) ± + EoffiO, s)] 


e 

Yh 


E Un,s - fN+l,s') 
Ar=l,s,s' 


F{N,s,s' 


H{N, s, s‘ 


rP 


x[S{hLO-EoH{N,s,s')) 
±5{nuj + EoH{N,s,s')] 


( 5 ) 


where the upper sign applies to Re(Jxx{^) and the lower 
sign to Imaxy{uj). It is clear from the form of equa¬ 
tion (5) that the optical matrix elements have restricted 
the possible transitions between LL to N^s N P s' 
with F{N,s,s') = {N,s\ja\N',s'){N',s'\jis\N,s)Tr and 
H{N,s,s') = -l + syi/4 + 2NP-syi/4 + 2{N + 1)P 
related to the energies of the possible optical transi¬ 
tions. While the complicated expressions for E{N, sp') 
specified in reference (26) (not repeated here as they 
are rather complicated and not particularly illuminat¬ 
ing) are to be used in the numerical results that we will 
present later, it is helpful to start with approximate ex¬ 
pressions which apply in the limit of small P. Defining 
E{N, 5, s') = E{N, 5, s')Eo/H{N, s, s') we find for the in¬ 
traband optical transitions to leading order (A^ = 1, 2,...) 
and for N = 


F{N,= -N/2{1-4:NP^)Eo (6) 


FiN, +, +) = -(1 + V)/2[l + 4(V + l)P^]Eo (7) 

with next corrections higher order in P and for the in¬ 
terband transition 

J^(V,+,-) = 2(l + 2V)p2^o (8) 


F{N, +) = -8V(1 + N)P^Eo (9) 

The first thing that needs to be emphasized is that the 
interband transition carry little optical spectral weight 
in the limit P 1. The (+, —) transitions go like the 
square of P and correspond to the absorption of pho¬ 
ton that have small energies as seen in Fig. 2 (short 
green arrows). Such transition vanish as P ^ 0. The 
(—, +) optical transitions (long green arrows) go like the 
cube of P and so are even less important. These cor¬ 
respond to finite photon energies with limiting value of 
2Eo as P ^ 0 as seen in Fig. 2. In contrast the intra¬ 
band transitions remain finite as P ^ 0 and so dominate 
the optical absorption. They represent the only possi¬ 
ble absorption processes in a pure non relativistic sys¬ 
tem with no spin-orbit coupling. This is clear from Fig. 
2 where the black arrows indicate the (+, +) intraband 
and red arrows the (—, —) intraband transitions. For a 
typical spintronic material at an applied magnetic field 
B of one tesla, Eq cy 1.16 meV for m = O.lrrie (with rrie 
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FIG. 3. (Color online) The dynamic longitudinal axx{(jo) 
(solid black curve) and transverse C7x2/(c^) (Hall,dotted red 
curve) conductivity in units of e‘^/h as a function of photon 
energy hw in meV. In all cases the Schrodinger mass m was 
set at its free electron value me, the magnetic field at H = 1 
Tesla, the broadening V — IK and the chemical potential fi 
is half way between the N — and iV = 1 landau level (LL). 
The Schrodinger magnetic energy scale Eq — 0.116 meV while 
the Dirac scale Ei is varied and is 0.116 meV for P = 1 in 
the top frame, 0.194 meV for P=2.8 in the middle frame and 
0.272 meV for P=5.5 in the bottom frame. The correspond¬ 
ing Dirac Fermi velocities are 4.8 x 10^ m/s, 8 x 10^ m/s and 
11.2 X 10^ m/s respectively. In the top frame arrows highlight 
the energies of some of the allowed optical transitions which 
contribute to the dynamic conductivity. The relevant LL en¬ 
ergy scheme is given in the inset where the chemical potential 
is shown as the horizontal dotted blue line. 


the bare electron mass) and for a Dirac Fermi velocity 
vp — 4.3 * lO^m/s, El 0.104 meV which corresponds 
to a P value of 0.008. Even if the magnetic field is 
reduced by a factor of fifty our expansion parameter P 
^0.4 and consequently the interband optical transitions 
are suppressed by a factor of (0.4)^ = 0.16 which is still 
small. It is important however to realize in this context 
that our expansion in small P implies that we are working 
at finite value of magnetic field (here 1 Tesla for definite¬ 
ness). Because Ei goes like and Eq goes instead like 
P, at very low field the Dirac energy scale will always be 
larger than the Schrodinger magnetic energy scale and 
we are mapped out of the small P regime. 


imation that P be very small and using the exact ex¬ 
pressions for the optical matrix elements in F{N, s, s') ^ 
In the top frame we show Reaxx{^) (solid black curve) 
and Im(Jxy{oo) (dotted red curve) in units of jh as a 
function of photon energy uj in meV. Here Eq = Ei = 
0.116meV, vp = 4.8 x 10^ m/s with chemical poten¬ 
tial /i set half way between the N = {) and N = 1 LL 
(/i = 0.17 meV). We have also used a residual scattering 
F = l/2r = IK which broadens out the LL although 
structures corresponding to the allowed optical transi¬ 
tions between levels are still seen and some of these are 
identified by arrows. Black and red are intraband, green 
interband in both main frame and on energy level inset. 
The inset gives the details of the allowed optical transi¬ 
tions including the position of the chemical potential at 
.17 meV (dotted blue horizontal line). We first note that 
En- is negative only for = 1 and 2 and is near the 
level with N = 0 for N = 3. It becomes larger than the 
value of the chemical potential for N = b which means 
that this level is unoccupied and hence the interband 
transition are limited to 7 with the —4 ^ 5 correspond¬ 
ing to a photon energy of .82 meV which is seen as a cut 
off for both longitudinal and transverse optical conduc¬ 
tivity as indicated by a green arrow. The lowest energy 
transition is intraband at 0.078 meV and the second in- 
traband is at 0.23 shown by the red and black arrows in 
the inset. Only 9 transitions are possible, two intraband 
and seven interband, one shown as the green arrow in 
the level scheme diagram. This is in sharp contrast to 
the pure Dirac case for which the interband transitions 
have no natural cut off in our continuum limit Hamilto¬ 
nian (1). Note that the magnetic energy scale is small 
and that, even for F = IK which is certainly reasonable, 
the broadening has eliminated the sharp peaks at the LL 
energies. Finally we note that, by choice, P was set to 
be 1 so that we are far away from the P ^ 1 limit con¬ 
sidered in the simplified Eq. (6) to (9). Nevertheless, the 
interband transition still correspond to suppressed opti¬ 
cal spectral weight. In the second frame of Eig. 3 we 
have increased the value of P to 2.8 and fi = 0.22 meV 
which is again half way between N = {) and N = 1 Lan¬ 
dau level. Now we see the emerging of a prominent and 
distinct interband background in Rea^xi^) on which is 
superimposed very broaden LL peaks. This background 
extends to almost ^2.0 meV and its height in units of 
/h is exactly tt/S (shown as the dotted blue line). Re¬ 
markably this is precisely the value found in graphene 
except for the factor of 4 which accounts for valley and 
spin degeneracy. The interband background is even more 
prominently developed in the lower frame of Eig. 3 for 
P = 5.5. Here the cutoff is < 3.0 meV. 


III. NUMERICAL RESULTS 

In Eig. 3 we show results of complete numerical cal¬ 
culations based on Eq. (5) without making the approx- 


We can gain insight into these results by comparing 
with results obtained in the topological insulator limit as 
we will do in the next section. Before we do this however 
it is also useful to consider the limit of zero magnetic 
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FIG. 4. (Color online) The pure Schrodinger quadratic in k 
dispersion curve (dashed line) compared with F^+(/c) (solid 
red) and E-{k) (solid blue) dispersion. The vertical arrows 
show the two extremum interband transitions that are possi¬ 
ble from which we identify kmin and kmax^ The shaded region 
emphasizes the area associated with such allowed transitions. 

field. In that case the formula for Reaxx{^) is given by 

- f{E+{k))]kdk 

p27T 

{hvFf d0TT5{hu) - E+{k) + E_{k)) (10) 

Jo 


Es(k) = - \-shvpk (11) 

2m 

The first thing to note about Eq. (10) is that the 
Schrodinger contribution to the energy drops out of the 
interband energy difference [E-^{k) — E-{k)] but it does 
remain in the thermal factor [f{E-{k)) — f{E-^{k))]. At 
zero temperature this factor reduces either to one in the 
interval for which interband transitions are possible or 
is otherwise zero. The kinematics involved are shown in 
Fig. 4. The red solid line gives E-\-{k) vs k, the blue 
is E_{k) and the black dashed curve, which is for com¬ 
parison, gives the dispersion curve when the spin-orbit 
term is zero i.e. vf = 0 in Eq. (11). The shaded re¬ 
gion shows the possible interband transitions and defines 
a minimum and a maximum value of momentum k^nin 
and /Cmax, which provide lower and upper limits on the 
integral in Eq. (10). These cut offs are due to the ther¬ 
mal factors which require an occupied initial state and an 
unoccupied final state. For graphene kmin corresponds to 
the onset of interband absorption at a photon energy of 
exactly twice chemical potential 2/i. In that case however 
kmax does not exist, since E-(k) is negative for all values 


of k and never crosses the chemical potential (/i > 0). 
Applying these restriction to Eq. (11) means that in¬ 
terband optical transitions are only possible for energies 
between 

^min — -^-|-(^min) (^min) 

and 

^max ~ -E-i-(A^max) E—(^kmax)' 

But the application of these cut off on momentum in Eq. 
(10) is the only effect of the thermal factors which, at zero 
temperature, have magnitude of one. This means that 
the remaining integral is identical to that for graphene 
leading to precisely the same absolute value of the uni¬ 
versal background. Of course now the background is re¬ 
stricted to the interval c^min to c^max- Thus Recjxxi^) for 
interband absorption is still equal to e^n/Sh but is non 
zero only for photon energies cJmin ^ ^ ^max- This 

remarkable result is consistent with the general trends 
found in our numerical work presented in Fig. 3 for the 
case of a finite B. Simple formulas for u;min and c^max are 
easily obtained 

Wmin = -2mv% + 2hVF\j^^Y + (14) 

Wmax = + 2hVF\J (^^^)2 + (15) 

We can check that in the limit m ^ oo (pure relativistic 
case) cjmax = oo and c^min = 2/i. In the opposite limit 
of Vf 0 both c^min and c^max are zero and there are 
no interband transitions. In all cases the width of the 
photon window over which the universal background has 
height e^ir/Sh is Co’max—^min = ^mv‘^p which is linear in m 
and quadratic in vf, and is independent of the chemical 
potential ji. However, both the onset and termination 
of the interband background do depend on the chemical 
potential /i which appears in the square root factor of 
Eq. (14) and this shifts both upper and lower limits on 
the interband absorption. A similar effect holds when a 
magnetic field is additionally applied, as we will describe 
below. First we note that the total optical spectral weight 
{Wib) in the interband background is Wjb = {e‘^7r/Shx 
Amvp). In addition there are the intraband transitions 
which give the Drude response. In the appendix we show 
that the optical spectral weight contained in the Drude 
when B = 0 is given by 

+ {2h^ / m)iJi — hvF)y/k?Vp + (2/i^/m)/i 
e^m{^lRvp + {2k? Im)[i + hvF)\JEJv\ + (2k? Im)[i 

where the first line comes from the (+, +) transitions 
and the second from the (—, —) transitions. Formula (16) 
applies also to the opposite limit of a topological insulator 
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FIG. 5. (Color online) Top frame, the dynamic longitudinal 
o'cccc(co’) (solid black curve) and transverse crccy(ct;) (Hall, dot¬ 
ted red curve) conductivity in units of e‘^/h as a function of 
photon energy hoj in meV. Parameters are the same as for the 
bottom frame of Fig.3 except that different values of chemi¬ 
cal potential are employed namely /x = 2.0 meV and jjl — 2. b 
meV. Bottom frame, the optical sum W(cUc) as a function 
of upper cut off cUc of the longitudinal conductivity axx((^)- 
The four curves are for fi — 0.28 meV (black), /x = 1.5 meV 
(red), 11 = 2.0 meV (blue) and /x = 2.5 meV (green). The 
dotted curves are for the case of H = 0 and are included for 
comparison with the solid curves for H = 1 Tesla. 


but in this case only the first line is retained since only the 
(+, +) transitions are possible. For the non-relativistic 
limit the Drude weight VFd is 

WD = ^[lJ.+ ]^mvl\. ( 17 ) 

Thus Wd does depend on the chemical potential /x in 
contrast to the interband contribution Wib- The total 
optical spectral weight Wt = Wd^Wib = 

We will return to this fact later when we consider spec¬ 
tral weight redistribution between intra- and inter-band 
transition as /x is varied and compare with the case of 
topological insulators. In the top frame of Fig. 5 we 
show results for the real part of the longitudinal dynamic 
conductivity {Re(7xx{^)i solid black curve) and the imag¬ 
inary part of the transverse Hall conductivity 
dotted red curve) for the same parameters as in the bot¬ 
tom frame of Fig. 3 but now the chemical potential has 
been increased. In the left frame /x = 2.0 meV and in the 
right frame /x = 2.5 meV. Increasing /x eliminates several 
interband transitions which now are no longer possible 
because of Pauli blocking. The optical spectral weight of 
these interband transitions has been transferred partly 


to the low energies intraband peak which becomes more 

prominent with increasing ji. In the lower frame of Fig. 

2 

3 a background of height ^ is well defined in the longi¬ 
tudinal conductivity down to photon energies < 1 meV. 
On the other hand in the top frame of Fig. 5 this back¬ 
ground is very suppressed in the region 1 meV to ~ /x 
where it is seen to jump to its universal value. For the 
case of 5 = 0 we saw that this threshold energy is c^min, 
defined in Eq. (14). This universal background then re¬ 
mains up to an upper cut off cjmax defined in Eq. (15) for 
5 = 0 at which point it drops to zero. Eor finite B these 
critical energies are some what different from those of Eq. 
(14) and vary with the magnitude of B. In Eig. 5 (top 
frame) B = 1 Tesla. Increasing B increases the distance 
in energy between the various Landau levels and so the 
energy of the minimum and maximum possible interband 
transition is changed. 

A more quantitative look at spectral weight distribu¬ 
tion between inter and intra band optical transition and 
its variation with value of the chemical potential /x is 
given in the lower frame of Eig. 5. There we present 
results for the total spectral weight contained in the lon¬ 
gitudinal conductivity below a variable upper photon en¬ 
ergy uj defined as 

nUJ 

W{(jj) = / Re(Jxx{^')duj' (18) 

Jo 

The units on W {uj) are ^ times meV’s and uj on the hor¬ 
izontal axis is in meV. Four values of chemical potential 
fi are shown /x = 0.28 meV (black), /x = 1.5 meV (red) 
/X = 2.0 meV (blue) and /x = 2.5 meV (green). Here 
the residual scattering rate F = liF, the Fermi velocity 
vf = 11200 m/s and the Schrodinger mass m is equal to 
the bare electron mass (rUe). 

The solid curves are for a magnetic field of B = 1 
Tesla while the color coded open circles which are shown 
for comparison, are for zero magnetic field. For the first 
black curve with small value of chemical potential /x, in¬ 
ter and intra band contributions to Reaxx{^) and to its 
integrated spectral weight W{uj) are not easily distin¬ 
guished although there is a clear onset, slightly below 
co’=0.5 meV, indicated by the black arrow, and there is 
also a cut off energy around 3.4 meV (see black arrow) 
beyond which the integrated spectral weight W(uj) ceases 
to increase. This is seen particularly clearly in the curve 
for 5 = 0 . For the other three cases considered, the 
intra and interband transitions are well separated. Be¬ 
low the onset of the interband transition at ujmin of Eq. 
(14), the B = 0 curve is completely fiat and contains the 
Drude spectral weight Wd of formula (17). At cj = uJmin 

there is an abrupt change in slope from zero to a finite 

2 

amount set by the value of the universal background . 
This kink, which we marked with a black arrow, while 
most pronounced in our B = 0 curves is also present at 
a slightly different energy in our B = 1 Tesla curves. 
This is expected and is traced to our numerical results 
for the absorptive part of the longitudinal conductivity 
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shown in the top frame of Fig. 5 where the onset of the 
universal interband background is clearly seen around 
uj = uJmin • The iutrabaud transitions below this on¬ 
set provide a peak in Re(Jxx{^) i displayed from cj = 0, 
which has however decayed to very small values before 
oj = ujjnin is reached. It is the displacement of the intra¬ 
band transition (which provides a Drude when B=0) to 
finite frequency in a nonzero magnetic field (B=l) which 
accounts for the gradual increase in W{uj) out of zero en¬ 
ergy {uj = 0) towards the Drude plateau of the B=0 case 
which is now never perfectly flat. 

Beyond the region of linear increase in IF(co’), in our 
B=0 case the optical spectral weight reaches saturation 
as can be seen in the black curves for which u;max of Eq. 
(15) falls within the range of energies shown in the figure. 
For the other cases one needs to go to energies beyond 
4 meV. The saturated value for B=0 is the sum of the 
Drude contribution of Eq. (17) plus the amount in the 
inter band which add up to + mvp]. 

This behavior is very different from what is known for 
the pure relativistic case such as for graphene for which 
no upper cutoff exists other than the Brillouin zone cut 
off. While we will see in a later section that modifications 
of the known graphene behavior also arise when a small 
subdominant Schrodinger piece is added to the dominant 
Dirac contribution, for the pure relativistic case spectral 
weight is simply redistributed between inter and intra 
contributions as /i is varied. Here we see that for materi¬ 
als with both Dirac and Schrodinger pieces in the Hamil¬ 
tonian there is no strict conservation of spectral weight 
with changes in /r, an issue we will return to in section 
V. While only /i changes in the various curves shown in 
the lower frame of Fig. 5, they do not merge as uj gets 
large. 


IV. COMPARISON WITH A TOPOLOGICAL 
INSULATOR 


scale given by the Dirac term only is Ei = 10 AmeV while 
the non relativity magnetic energy scale for m = rrie is 
0.116 meV, two orders of magnitude smaller. Secondly 
the interband optical peaks in Reaxx{^) remain to very 
high energies as in the case of graphene previously dis¬ 
cussed in Ref. (17). Just as we found it enlightening in 
Fig. 3 to consider several values of vp for fixed value of 
m, here we fix and consider 3 values of m. For m = rrie 
third frame from top in Fig. 6 the results are not signif¬ 
icantly different from those in pure graphene. As m is 
decreased to O.lme a splitting of the peaks in Re(Jxx{^) 
into pairs is seen. This arises because the Schrodinger 
term breaks the particle hole symmetry of the pure rel¬ 
ativistic case and the optical transition —N to N 1 
no longer has the same energy as —{N + 1) to V. This 
splitting is even more pronounced when m is decreased 
further to m = 0.05me as in the top frame. In this case 
not all peaks can be easily identified as split pairs. Of 
course even for m = rUe , (third frame down) there should 
in principle be a splitting but here we have taken a smear¬ 
ing parameter F = 15A which is enough to merge them 
so that a single peak is effectively seen. 

We will return to the results of Fig. 6 later. First it 
is useful to consider analytic results which we can obtain 
only in the limit of P~^ 1 when the non relativistic 

correction to pure Dirac is small. In this case the opti¬ 
cal spectral weight in units of e‘^/2h which is associated 
with the allowed transitions at zero temperature have the 
approximate for m^^i^^ 


w-,-) = [ 


vW+VVTT 


4^2 


1/Vp 


16 [v/iV(l+iV)- 2 ^^^g^ 


-]Ei 


(19) 


which will not be needed for /i ^ 0 because such transi¬ 
tions are Pauli blocked 


We next compare the results of Fig. 3 with results 
obtained in the topological insulator limit shown in Fig. 
6. The three top frames give the absorptive part of the 
diagonal conductivity Reaxx{^) (solid black curve) and 
comparison with the imaginary part of transverse con¬ 
ductivity Im(Jxy{uj) (dotted red curve) in units of e^/h 
as a function of photon energy uj in meV. We have taken 
the Fermi velocity to be 4.3 x 10^ m/s, the magnetic 
field to be 1 Tesla, the scattering rate F = 15A and the 
Schrodinger mass equal to the bare electron mass {rrie) in 
the third lowest frame, m = O.lme in the second lowest 
frame and m = 0.05me in the top frame. Decreasing m 
moves us further away from the pure Dirac case which 
would correspond to graphene. In all these frames the 
chemical potential was set to fall halfway between N = 1 
and N = 2 Landau level. The first thing to note about 
these results is that the magnetic energy scale associated 
with the Landau levels is much larger than that in Fig. 3 
for the non-relativistic limit. Here the relativistic energy 


T{N,+,+) = -[ 


4x/2[-W +/VTT] 
^E 

i+2Ar 


'l6WN{l + N)-2^q^A' 


E{N,-,+) = [ 


-1 


4V2 [Vn + vrrTi] 
iIVp 




]Ei 


T{N,+,-) = 


4V2[VN + y/WTl] 
iIVp 


16 [^Nil+N)+ 2 Jq^ 


:]Ei 


( 20 ) 


( 21 ) 


( 22 ) 
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where we are working to lowest order in \j\fP. This 
provides the first correction to the pure relativistic limit. 
Note that, in addition to the T{N,s,s')^ the absorptive 
part of the longitudinal optical conductivity Re(Txx{^') 
depends additionally on thermal factors which give a 
+1 or —1 atT = 0. consequently the optical spectral 
weight under an allowed optical transition between LL 
which must be positive, is in all cases equal to the ab¬ 
solute value of J-{N^ s'). This means that the spectral 

weight corresponding to ( 20 ) and ( 21 ) carries an addi¬ 
tional minus sign so that the optical spectral weight asso¬ 
ciated with ( 21 ) is reduced over its pure Dirac limit value 
while that associated with ( 22 ) is increased by exactly 
the same amount. When this correction is dropped, we 
recover the result of Gusynin, Sharapov and Carbotte^^ 
for graphene except for a missing factor of four account¬ 
ing for spin and valley degeneracy not present for topo¬ 
logical insulators. Taking the limit of large N the op¬ 
tical spectral weight contained in the intraband tran¬ 
sitions is given by|^ J’(A^,-h,+)Ar^iarge- Noting that 

-— j=^^—== 2\fN and that — , ^ — 2 ^ 

[-vW+vW+T] [7Ar(i+Ar)-2 ^/^,y^ ] 

we get 

“ +)Ar^large = -h l/\/P 

But + ^o(^)^ ^(1 - ^) 

so that 

+ ( 24 ) 

2 

In the limit of m = oc this reduces to which is ex¬ 
actly the amount of optical spectral weight there would 
be under the Drude in the 5 = 0 limit in graphene except 
for the factor of 4 in the denominator which would be can¬ 
celed by a degeneracy factor for two spins and two valleys. 
When m is large but not infinite the spectral weight un¬ 
der the intraband transition line is increased from a nor¬ 
malized value of 1 by an amount ^^^2 • For our expan¬ 
sion to be valid we still need . ^ 2 to be small which is a 

A.mVp 

more restrictive condition than simply \/\fP ^ 1. Nev¬ 
ertheless it shows clearly that the optical spectral weight 
residing in the various optical lines is changed from the 
pure Dirac case when a subdominant Schrodinger piece 
is also present in the Hamiltonian. 

In the large N limit the intraband line is the cyclotron 
resonance line of semiclassical theory. The cyclotron fre¬ 
quency was worked out in Ref. [26] and found to be to 
lowest order correction for a small Schrodinger contribu¬ 
tion. 


hiOc = 5[1 + ^Ij/{mvjr)] (25) 

So that hwc is increased for m 7 ^ oc as is the spectral 
weight under this line. It is interesting to compare the 


spectral weight of the cyclotron resonance line with the 
Drude weight (Wd) for the zero magnetic field case. The 
expression for Wd is (Eq. (16), first line only) 

e^m[h?v% + {2Y? jm)^ — hvF\/^v\ + {2}^ jm)^ 

8^3 ( 26 ) 

which is valid for any value of m. Assuming m to be 
large but not infinite we can expand (26) and obtain a 
first correction to pure Dirac, we get 


Wd 


e^TT 


fi[i 


2mvp 


(27) 


which agrees perfectly with the spectral weight under the 
semiclassical cyclotron line in this approximation. 

In the comparisons made so far between the pure Dirac 
case and a topological insulator with finite Schrodinger 
contribution we have considered the chemical potential 
as fixed. It is important to realize that this does not 
correspond to a fix doping as m is varied. For a fix density 
of charge carriers (n) away from the neutrality point we 
have that 



I —mvF 
An h 


mvF 

h 



(28) 


where we have used kmin determined from Fig. 4. Eq. 
(28) holds whatever may be the value of m. For large m 
retaining the first leading correction we get 


Af hVF 


mvp 


(29) 


from which it follows that 


/i AFnhvF{'^ + 


mvF 


(30) 


For m = oc we recover the known result for pure Dirac 
fermions ji = ^AfuHvf- For m not infinite, /i is in¬ 
creased over the pure relativistic limit value because of 
the change in the electronic dispersion curves. This in¬ 
crease in chemical potential is a direct consequence of 
the narrowing of the cross-section of the conduction band 
cone due to the subdominant Schrodinger term. 

Returning to our approximate equations for 
5(y+,+), 5(A^,+,—) and E{N,—,-\-) we note 

that in all cases the corrections for a non zero value of 
1/P is of order l/\fP and more importantly, that the 
spectral weight associated with the interband transi¬ 
tion remains finite even when I/a/P = 0. This is in 
sharp contrast to what we found in the non-relativistic 
limit where they vanish when no subdominant Dirac 
correction is included. For the pure relativistic case 
the photon energy associated with both (+, —) and 
(—,-h) interband transitions is the same. But here with 
\/\fP 7 ^ 0 the (+, —) transition is slightly shorter than 
is the (—,+) transition. We also see, noting the sign of 
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the thermal factors in Eq. (5), that the optical spectral 
weight associated with the larger photon energy has 
the smallest spectral weight. To first leading order in 
1/a/P however the total spectral weight under the two 
split lines is unchanged from the pure relativistic case 
where \/\fP = 0. In graphene it is well knownii that 
as /i is increased so as to cross the N\h. LL line only 
one half the spectral weight remains in the N^th line 
and all others for n < have disappeared with the 
entire spectral weight lost in the interband transition 
reappearing in the single intraband line which moves 
to lower photon energy with increasing fi and picks up 
more intensity. This sum rule on the redistribution 
of the spectral weight between inter and intra lines is 
encapsuled in the equationii 

^ V^rTT+vWTT +Vif ~ VFTT- 

(31) 

derived in Ref [17], see their Eq. (33). The first term 
on the left is related to the spectral weight associated 
with all lines that have completely disappeared while the 
second is half of the spectral weight of the line n = N. 
The right hand side is the spectral weight in the intra¬ 
band line. It is clear from our approximate Eq. (19-22) 
and Eq. (5) that this sum rule will no longer hold when 
1/Vp^o since pairs of lines that have completely dis¬ 
appeared will have the same spectral weight as for the 
pure Dirac case while the one for n = N will have a cor¬ 
rection of order 1 / ^/P (reduction) which does not match 
precisely the intraband correction of order 1 /\/P. 

Next we return to our numerical results of Eig. 6 which 
do not require the \/\fP 0 limit and these show ex¬ 

plicitly a change in optical sum with Schrodinger admix¬ 
ture. The first peak in each of the three top curves for 
Reaxx{uj) come from the intraband optical transitions 
while all other peaks are interband. As the Schrodinger 
contribution is increased through a decrease in m all in¬ 
traband lines are seen to shift slightly to higher energies 
as compared to the pure Dirac case (which close to the 
m = rUe curves). Also the optical spectral weight con¬ 
tained in the intraband line increases as the Schrodinger 
mass is decreased. This fact is more easily seen in the 
lower frame of Eig. 6 where we give results for W{uj) of 

Eq. (18) as a function of uj in meV with W in units of 

2 

^meV. All the curves start from zero and rise rapidly as 
we integrate over the intraband line after which it has a 
plateau which is nearly but not quite constant, followed 
by a more rapid rise modulated by small wiggles that 
reflect the LL structure. These wiggles would of course 
be more pronounced if we had reduced the broadening 
r. Here it is 15i^. The height of the first plateau in 
W{uj) gives the amount of spectral weight contained in 
the intraband transition and this clearly increase as we 
move further away from the pure relativistic case. Differ¬ 
ences in W{uj) remain to high energies and these reflect 
the small Schrodinger admixture which we have added to 
the Hamiltonian. 



FIG. 6. (Color online) Top frame is the dynamic longitudinal 
(solid black) and transverse axy{uj) (Hall, dotted red) 
conductivity in units of e^//i as a function of photon energy huj 
in meV. The Dirac fermi velocity is set at i’f = 4.3 x 10^ m/s 
which is nearly two orders of magnitude larger than used in 
the top frame of Fig. 3 and is representative of topological in¬ 
sulators. Top, second and third frames are for the Schrodinger 
mass equal to me, O.lme and 0.05me respectively with rUe 
the bare electron mass. The bottom frame gives the optical 
spectral weight W(ijJc) in units of e^/hxmeY as a function of 
photon energy hw in meV. The curves are for H = 1 Tesla. 
The solid black is for m = me, dotted blue for m = O.lme 
and dotted red is for m = 0.05me. 


V. REDISTRIBUTION OF SPECTRAL 
WEIGHT WITH /x 

Next we return to the issue of how optical spectral 
weight gets redistributed between inter and intra band 
when the chemical potential fi is varied i.e. a change in 
doping, without assuming 1 /a/P << 1. We begin with 
the case of H = 0 because in that instance, as we have 
seen, we can get analytic results. In previous sections 
we showed that the Drude weight is given by Eq. (16). 
Only the first line in this equation applies in the rela¬ 
tivistic limit while both lines (Eq. (16)) contribute in 

the non-relativistic case. Eurther, we saw that the in- 

2 

terband transition provide a constant background of 
in magnitude, with sharp absorption edge at c^min which 
terminates at cjmax (see Eq. (14) and (15)). The total 
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0 2 

optical spectral weight in this background is Amvp^. 
For comparison with the relativistic limit let us begin 
by computing the optical spectral weight that is miss¬ 
ing from the universal background between cj = 0 and 
uj = cjinin- For m ^ oo we can expand the expression for 
c^min and obtain to lowest order 


which gives the known result that the interband absorp¬ 
tion edge starts precisely at 2 /i in the pure relativistic 
limit. When m is not infinite the absorption edge has 
shifted to an energy somewhat less than 2/i. The missing 
optical spectral weight is 




-JL.] 

2mvp 


(33) 


2 

which is also somewhat less than /i found in the limit 
m = oo. Now in the limit of a topological insulator the 
Drude weight is given to leading order in an expansion 
for large m in Eq. (27) and is not equal to the missing 
weight in the universal background given in Eq. (33). 
So the spectral weight redistribution sum rule which is 
operative in the pure relativistic limit breaks down when 
a subdominant Schrodinger piece is added to the Hamil¬ 
tonian. 

The opposite limit (non-relativistic case) is also of in¬ 
terest. In that case both terms in the Eq. (16) need 
to be retained and to lowest order in vf^ Wd is given 
by Eq. (17). Eurther the missing weight in the universal 
background below the interband absorption edge is in the 
same approximation 


(34) 

which does not depend on /r in the same way as (17). It 
is clear that no sum rule applies in this limit as well. The 
pure relativistic case is unique and special. 

Similar results hold when a magnetic field is present as 
we can see in Eig. 7 which gives the total spectral weight 
below UJ (IF(co’)) under the real part of the longitudi¬ 
nal conductivity Reaxxi^) in units of (e^/h)xmeV as a 
function of uj in meV. Eour cases are shown with B = 1 
Tesla, F = 15i^, vp = 4.3 x lO^m/s and the Schrodinger 
mass m = O.OSmg. The four values of chemical poten¬ 
tial are (solid black curve) /i = 29.1meV, (dashed green) 
fi = 35.7meV, (dashed dotted red) fi = 41.8meV and 
(blue dots) = 47.4meV. For the first curve the open 
circles compare with the B = 0 case which starts at a 
finite value at cj = 0 because we have assumed the clean 
limit and so the intraband transitions contribute a Dirac 
delta function to the frequency dependent conductivity. 
W{uj) then remains completely flat until the onset of the 
interband transitions set in at c^min = 47.8meV given by 
Eq. (14). This is quite a lot smaller than the value of 
2fi 60meV. The curve for B = 0 follows well the solid 





b 

0 


O) 

0 





CO (meV) 


FIG. 7. (Color online) The partial optical spectral weight 
W (cj) to an upper cut off uj (in meV) for four values of chem¬ 
ical potential /x = 29.1 meV (solid black curve), /x = 35.7 meV 
(dashed green curve), /x = 41.8 meV (dash-dotted red curve) 
and /X = 47.4 meV (dotted blue curve). In all cases the mag¬ 
netic field is B=1 Tesla, the Dirac fermi velocity is 4.3 x 10^ 
m/s and the Schrodinger mass is 0.05?7Xe. The open circles are 
for comparison with the solid black curve but now B is set to 
zero. The light dashed line extrapolates the straight line de¬ 
fined above the interband threshold in the solid black curve. 
Its intercept with the vertical axis is a measure of the devia¬ 
tion from pure Dirac behavior brought about by the presence 
of a subdominant Schrodinger piece in the Hamiltonian (1). 


black curve for H = IT except of course that in this in¬ 
stance the Drude is replaced by the intraband peak and 
so W{uj) starts from zero at cj = 0 and becomes equal 
to the black dots only after all the spectral weight in 
this line has been picked up. Above the interband on¬ 
set both curves rise linearly with increasing uj with the 
slope of this line set by the value of the universal back¬ 
ground. If we extrapolate it to u; = 0 we find that it 
intercepts the uj = 0 axis at a non zero value given by 

^ujn,in{A - 1) = ^mvpA - If with A = ^1 + ^. 

For m = oc, A = 1 and this intercept is equal to zero 
which applies only for the pure relativistic case. By con¬ 
trast in topological insulators A is not one. Should one 
know independently the value of the chemical potential 
one could use the value of the intercept to 

determine mv^ which is a central parameter in assessing 
the deviations from the pure relativistic case expected in 
a given topological insulator. Note with reference to the 
lower frame in Fig. 6 that an extrapolation of the near 
linear region of the black curve for m = rUe, to cj = 0 will 
give almost zero intercept and this reflect the fact that 
in this case the Schrodinger correction to pure Dirac is 
small. 

Returning to the data in Fig. 7, W{uj) rises to a larger 
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magnitude for the plateau betweeu the saturatiou of the 
iutrabaud spectral weight aud the threshold for the iu- 
terbaud trausitiou as the chemical poteutial is iucreased. 
This reflects the fact that as /i iucreases, more iuter- 
baud hues are lost because of Pauli blockiug aud some of 
their optical weight is trausferred to the iutrabaud hue. 
The ouset of the iuterbaud trausitiou beyoud the iutra¬ 
baud plateau iu W{uj) iucreases with iucreasiug /i but 
the liuear rise above this critical photou euergy has the 
same slope iu all cases. The small variatious about a 
straight hue seeu iu this regiou reflect of course the uu- 
derlyiug Laudau level structure which is smeared almost 
completely by the broadeuiug T = IbK. At higher eu- 
ergies all the curves become parallel to each other but 
are displaced upward as /i is iucreased. Iu grapheue, as 
cau be seeu iu Fig. 7 of Ref. [17], the curves would all 
merge iu this regiou iu coutrast to what we hud here for 
topological iusulators wheu the value of the subdomiuaut 
Schrodiuger part to the Hamiltouiau is siguiflcaut. It is 
the value of the iutercept with the verticle axis of the 
straight hue variatiou at large euergies which sets the 
displacemeut betweeu these parallel hues aud as we have 
seeu grapheue correspouds to the limit of zero iutercept. 

VI. SUMMARY AND CONCLUSIONS 

We have calculated the absorptive part of both lougi- 
tudiual aud trausverse (Hall) AC dyuamic magueto cou- 
ductivity for quadratic electrouic bauds modified by a 
small spiu orbit coupliug. This provides a mechauism 
for iuterbaud optical trausitious betweeu the split heli¬ 
cal bauds. These optical trausitious are additioual to the 
maiu iutrabaud (Drude) absorptiou which is the ouly pro¬ 
cess possible wheu the spiu-orbit coupliug is uot preseut. 
We hud that for the rather small values of Fermi veloc¬ 
ity represeutative of the semicouductors preseutly used 
iu spiutrouic applicatious, the iuterbaud absorptiou is 
small aud that this backgrouud is likely uot to be seeu 
as distiuct from the maiu Drude iutrabaud absorptiou. 
However as the Fermi velocity vf is iucreased a distiuct 
iuterbaud backgrouud emerges which takes ou the same 
uuiversal values of ^ as observed iu graheue except for 
a degeueracy factor of 4 uot applicable iu this case. The 
uuiversal backgrouud however is couflued to a defluite 
very limited rauge of photou euergies with upper aud 
lower cut off depeudeut ou the value of the chemical po¬ 
teutial aud also ou the maguitude of the maguetic field 
wheu B is also preseut. Comparisou of these results with 


those obtaiued iu the opposite limit, wheu iustead the 
Schrodiuger term is a small perturbatiou ou a domiuaut 
spiu orbit term, provides uew iusight iuto the case of 
topological iusulators. There, the iuterbaud trausitious 
are much more promiueut, exteud to large euergies aud 
are cut off ouly by the baud edge. For uo maguetic field 
{B = 0) we were able to obtaiu aualytic results which 
are valid for auy value of m aud vf aud these coufirmed 
aspects of our fiuite B results. We showed that the iuter¬ 
baud backgrouud coutaius a total optical spectral weight 
of which is iudepeudeut of ji. This spectral 

weight is small wheu either m or vf is small aud it is 
large wheu both m aud vf are large as iu topological iu¬ 
sulators. Iu particular for m ^ oc we recover the kuowu 
result for grapheue where ouly the Dirac coutributiou is 
preseut, aud the iuterbaud trausitious exteud to the baud 
cut off. 

Wheu the chemical poteutial /i is iucreased, a well 
kuowu result for grapheue is that all the optical spectral 
weight lost iu the iuterbaud trausitious is fouud traus¬ 
ferred to the iutrabaud trausitious. This couservatiou of 
spectral weight uo louger applies wheu a subdomiuaut 
Schrodiuger piece is added to a pure Dirac Hamiltouiau. 
While these violatious are small for typical parameters 
associated with preseut day kuowu topological iusula¬ 
tors, measuriug them allows oue to determiue the crit¬ 
ical parameter mv‘^ which provides the iuformatiou ou 
the expected deviatious of optical properties from the 
pure Dirac case. Iu the spiutrouic limit we hud that the 
amouut of spectral weight iu the iuterbaud backgrouud 
is almost fixed while at the same time the iutrabaud cou¬ 
tributiou iucrease with iucreasiug chemical poteutial. 
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Appendix A: Derivation of the Drude spectral 
weight at zero magnetic field 

The Hamiltonian is given by Eq. (1) and here we write 
a = hvF‘ 


2m 

The velocity is 

dHo n^k 


Hq ^ T (^{kxO’y kyCJx) 


Vr = 


dkx m 
dHo h^k 

dky 


m 


Jy ruy^ 

+ OLCTy = V]jl + dCfy^ 

^ — acjx = vdI — OiCFx' 


(Al) 


The matrix Green’s function can be expanded in matrix 
spectral function A(k, uj) as 


Gik,z) = j' 


duo A(k, uj) 
27r z — UJ 


(A 2 ) 


The conductivity crxx{(^) is given by 

„2 .oo ,oo [f{co0-f{0J2)] 

,_oo j-oo 2Tr UJ - UJ2 + UJI + i5 

V Tr{vxA{k,uJi)vxA{k,uj2)) (A3) 


tuj 2 -k J_, 


and the real part of the conductivity at zero temperature 
is given by 


Re(Txx{i^) = 


2uj dTT^ 


f'k'cut 

kdkdO 


r duji 


Tr{vxA{]^,uji)VxA{]^,uj + cji)) (A4) 

The spectral function is the imaginary part of the Green’s 
function 

A(k, x) = — 2 /mG(k, x) (A5) 

The Green’s function is given by 

Go{k,iujn) = :^ ' <^)Go(k, s, iw„) (A6) 


where 


and 




Ffc = 


( hy , kx 1 0) 


Go(k,s,iw„) = 7 


+ yu - - sak 


(A7) 


(AS) 
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After taking the trace we get 


duji 


Reaxx{uj) 

2 r^cut 

— / kdk 

4^ Jo j-bj 

27ra‘^[6{uji — ak)5{uj2 — oik) + 5{uji + ak)5{uj2 + oik) 

+ ak)5{uj2 — Oik) + S{uJi — ak)S{uj2 + Oik)] 

+27r—^-[^(co’i — ak)5{uj2 — Oik) + ^(cji + ak)5{uj2 + Oik)] 

/^2 _ _ _ _ 

-\-4:7ra -[^(cJi — ak)S{uj2 — Oik) — S{uJi + ak)S{uj2 + Oik)] 

m 


|/'(fco)| = I - + q;| = a/q !2 + 2 h^/m{u)i + /x). Thus 

'kcut 


rf^cut 

/ dk • k 6 {uji + ak) 5 {uj 2 + oik) 

Jo 

m a — o? + 21 i?/m{uJi + /i) 


k? \/oL^ 2 h ?+ /i) 


5{uj) (A14) 


where 


CJi = CJi + /i - 


kJk‘^ 


2 m 


6^2 = CJi + CJ + /i - 


li?k‘^ 

2 m 


(A9) 

(AlO) 


/ dk • k‘^ 6 {uji + ak) 6 {uj 2 + c^k) 

Jo 

= + ,415) 

d JcP- + 2h?/m{u!i + n) 

pl^cut 

/ dk • k^ 6 {uji + ak)5{uj2 + Oik) 

Jo 

= + ,415) 

k o? + 2h?lm{uJi + /i) 


and hence 


The term with S{uJi-\-ak)S{uj 2 — Oik) and S{uJi — ak)S{uj 2 -\- 
ak) will be zero for topological insulator but nonzero 
for spintronics. For topological insulators we have for 
/i > 0 ,cl;i ^ 0 ,cji + /i — — ak = 0 ^ k = ■p-[— 

^ya‘^ + 2 h?/m{u)i + /i)], and 

r^cut 


f 

J —u 


duJi 


a — ya^ + 2^?777i(c(Ji+70 


-w + 2 h?lm{uJi + /i) 

acj 

— CJ 


+ { 2 lT?/m)ii 


Similarly 


pi^cut 

/ dk • — ak) 6 {uj 2 — oik) 

Jo 


m —a + + 21 i?/m{uJi + /i) 

+ 2 Ji?/m{uJi + /i) 


.0 

/ dcji 

J —UJ 


S(cu) (All) 


(a — o? + 2k? Im{uj\ + /i))^ 
+ 2 h‘^/m{uoi + /i) 

(a — + { 2 h? f m)UJ 


f 


and 


dk • /c^^(co’i — ak)S{uj 2 — ak) 


= ," ,3 (-« + + 2BVm(3,i +1.))^ ,, ^ 

k 2 h?/m{uoi + jn) 


L 


° , {a - Jo? -V 2kfl/m{oJi + n))^ 

dull- 

— UJ 


^ya‘^ 2ffi/m{uoi + fi) 




^/a‘^ ^ { 2 Ji?/m)iJ 


Thus for /i < 0, 


e‘^m{a - ^ya‘^ ^ ( 2 k?lm)[})^J + ( 2 k?lm)[x 
Reaxx[^) = - ^ 


and for /i > 0 , 


(A17) 


rkcut 

/ dk • k^ 6 {uji — ak)5{uj2 — Oik) 

Jo 

= '3 + 

k^ + 2 JR /m{uji + /i) 

For fi < 0,co’i 0,co’i + /r — ^—h = 0 ^ /co = 

t[a - + 2fi»/m(3., + rt]^ Because t„.. = ^ and + (2hym)^ - n)+ (2llV^5 

the energy at this point F’min = ~^2 • The delta func- (A18) 

tion can be rewritten as S{uJi + ak) = where For spintronics we can perform similar algebra. 



